Edge states and the integer quantum Hall conductance in spin-chiral ferromagnetic 

kagome lattice 
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We investigate the chiral edge states in the two-dimensional ferromagntic kagome lattice with 
spin anisotropies included. The system is periodic in the x direction but has two edges in the y 
direction. The Harper equation for solving the energies of edge states is derived. We find that 
there are two edge states in each bulk energy gap, corresponding to two zero points of the Bloch 
function on the complex-energy Riemann surface (RS). The edge-state energy loops parametrized 
by the momentum k x cross the holes of the RS. When the Fermi energy lies in the bulk energy gap, 
the quantized Hall conductance is given by the winding number of the edge states across the holes, 
which reads as cr£y SC = — ^-sgn(sin 4>), where <j> is the spin chiral parameter (see text). This result 
keeps consistent with that based on the topological bulk theory. 

PACS numbers: 73.43.-f, 73.43.Cd, 71.27.+a 
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Recently the quantum transport of electrons in spin- 
orbit coupled 13, 0, Q or spin-chiral ferromagnetic sys- 
tems 0, la, 0, 0| has been a focus of intense interest in 
condensed matter physics. One typical spin-chiral ferro- 
magnetic system is represented by pyrochlore compounds 
-R2M02O7 (i?=Nd, Sm, Gd), in which the spin configu- 
ration is noncoplanar and the spin chirality appears. As 
a consequence, the quantum transport of electrons, espe- 
cially the transverse conductivity a xy , is expected to be 
affected by the presence of spin chirality. Ohgushi et al. 
Q have first pointed out that the chiral spin state can 
be realized by the introduction of spin anisotropy in an 
ordered spin system on the two-dimensional (2D) kagome 
lattice, which is the cross section of the pyrochlore lat- 
tice perpendicular to the (1,1,1) direction |9(. In this 
case, it has been shown in the topological bulk theory 
[1, Hfll that the presence of chiral spin state may induce 
gauge-invariant nonzero Chern number, thus resulting in 
a quantized Hall effect in insulating state. 

In this paper we turn to study the 2D kagome lattice 
with two edges, which, as will be shown below, displays 
two chiral (instead of nonchiral) edge states localized 
near the sample boundaries. Closely followingthe topo- 
logical edge theory established in last decade [IJ, Qjl, Hal > 
we first derive the transfer matrix (namely, the Harper 
equation [2 HE, [lfj ) for solving the energies of the edge 
states. Although the transfer matrix elements are no 
longer polynomials of the energy e for a fixed spin chiral 
parameter 4> (except for special cases 0=0, ±7r/2, 7r), we 
find that the transfer matrix method is also applicable 
in this system. Then by numerical calculation, we find 
that there are two edge states in each bulk energy gap, 
corresponding to two zero points of the Bloch function on 
the complex-energy Riemann surface (RS). Remarkably 
different from the case of a square lattice under a mag- 
netic field [ll|, the edge-state energy loop moves across 
(not around) the holes in the RS in the present model. 
The two edge-state energy loops lying in the same energy 
gap are tangent at one point and their appearance shows 
a "00" structure. Furthermore, we obtain that when the 



Fermi energy lies in the bulk energy gap, the quantum 
Hall conductance given by the winding number of the 
edge state can be written as er°y SG =— ^-sgn(sin</>). This 
result based on the topological edge theory keeps con- 
sistant with that based on the topological bulk theory 

Following Ref. Q, we consider the double-exchange 
ferromagnet on the kagome lattice schematically shown 
in the left panel in Fig. 1 [U, Here the triangle 

is the one face of the tetrahedron, and the easy axis of 
the spin anisotropy points to the center of each tetrahe- 
dron and has an out-of- plane component. In this situ- 
ation the three local spins on sites A, B, and C in the 
left panel in Fig. 1 have different directions and the 
spin chirality emerges. The effective Hamiltonian for the 
hopping electrons strongly Hund-coupled to these local- 
ized spins is given by H — J2nn ^ij c i c i w ^h — 
t(xi\Xj) — te laij cos -4f. Here the spin wave function \xi) 

is explicitly given by \xi) — [cos^, e % '^ i sin^] T , where 
the polar coordinates are pinned by the local spins, i.e., 
(Xi\Si\Xi) = \ (sintfjcos^,, sin sin 4, cosi?*). % is 
the angle between the two spins Si and S^. The phase 
factor a,j can be regarded as the gauge vector poten- 
tial ci/i(r), and the corresponding gauge flux is related to 
scalar spin chirality Xijk=Si-(SjXSk) I n periodic 

crystal lattices, the nonvanishing of the gauge flux re- 
lies on the multiband structure with each band being 
characterized by a Chern number [l8|, Following 
Ref. [1, [1(3] we set the flux originated from the spin 
chirality per triangle (see Fig. 1) as <f>, which satis- 
fies = e l ( a AB+a B c+a C A) _ xhe flux penetrating one 
hexagon is determined as —2c/). We take the gauge, in 
which the phase of i^j; is the same for all the nearest- 
neighbor pairs with the direction shown by the arrows in 
the left panel of Fig. 1. It should be pointed out that 
the net flux through a unit cell vanishes due to the can- 
celation of the contribution of the two triangles and a 
hexagon. Also noted is that the time-reversal symmetry 
is broken except for cases of 4>—0,tt. In the following we 



2 




This is the Harper equation [IJ, 0j| . Equation ([3]) can 
be represented in the following matrix form: 



where 



V(m+1)A \ _ jCf ( f ) ( ^mA 
^mA ) ( ' \ *{m-l)A 



M(e) 



P -1 

1 



(4) 



(5) 



FIG. 1: (Color online) Left panel: Two dimensional spin- 
chiral ferromagnetic kagome lattice. The dashed line repre- 
sents the Wigner-Seitz unit cell, which contains three inde- 
pendent sites (A, B, C). It is assumed that each site has a 
different spin anisotropy axis. The arrows means the sign of 
the phase of the transfer integral tij. Right panel: The 2D 
kagome lattice system with edges along the y direction. 



and P= 2cos(fcl/2)( e 4 + ) 2cos0) " 2C0S (¥)■ In the following 

we do not explicitly write the subscript A in Eq. f4|. 
Then we get a reduced transfer matrix linking the two 
edges as follows 



M(e) 



(G) 



change notation i — ► (Ims), where (7m) label the kagome 
unit cell and s denote the sites A, B and C in this cell. 
The size of the unit cell is set to be unity throughout this 
paper. 

Now we suppose that the system is periodic in the x 
direction but has two edges in the y direction (see the 
right panel of Fig. 1). Since the system is periodic in 
the x direction, we can use a momentum representation 
of the electron operator 

1 Ax {lms)lms{kx) ^ (1) 



-•(Zms) 



where (Xn ms \,Ya ms \) are the coordinate of the site s 
in the unit cell (Im) and k x is the momentum along 
the x direction. Let us consider the one-particle state 
l*(M)=Em S ^ms(k x )^ ms (k x )\0). Inserting it into the 
Schrodinger equation H\^)—e\^) : we can easily obtain 
the following three eigenvalue equations for sites A, B, 
and C, 



+ e'3 



(m+l)C 



e$ mB = 2e 4 3 cos ( — ) * mC 



+ e'3 



e ^ k ^ m A + e l ^ k ^ (m _ 1)A 



k. 



(2) 



Eliminating the B- and C-sublattice sites, we obtain the 
difference equation, 



, ! 1, f 1 • cos 2 ' / ''' 



8cos 2 (y ) cos0 



mA 



= 2 < | y ) (e + 2 cos <jy) [^ (m +i)A + ®( m -i)A\ ■ (3) 



where 



M(e) = 



M(e) L «- f Mll(e) M »(^ (7) 
M ^ " {M 21 (e) M 22 (e) J' Uj 



For general 0, which varies in a range between — 7r and n, 
Mij (e) (i, 7=1,2) are not polynomials of e. At four spe- 
cial values, i.e., 4>=0, ±-|, tt, however, they can be writ- 
ten as polynomials of e, with the degree of 2L y for Mu, 
2L y — 1 for Mi 2 and M 2 i, and 2-Lj, — 2 for M 22 . In fact, 
in the spin-chiral cases of <^=±^, the factor p in M^ is 

reduced to P= 2co £ s (fc 4 / 2 ) — 2 cos (4^). Whereas, in the 
spin-nonchiral cases of 0=0, n, the factor p in M^ is re- 



2 cos I 



All kinds of solutions 



duced to p=- 2cos(fca/2) ~ ~« v „ 
from Eq. ([6]) are obtained by different choices of and 

Now we investigate the energy spectrum of the one- 
dimensional problem with special attention to the edge 
states. The boundary condition of this problem is 



(8) 



With Eqs. (jGj) and (0, one can easily obtain that the 
solutions satisfy 



Mai (e) = 0. 



(9) 



From 



= M 2 (e) 



*L„-2 



= M(e) 







one can find that 



-Mn (e)*i 



(10) 



If we use a usual normalized wave function, the state is 
localized at the edges as 

|Mn (e) | <C 1, localized at y « 1 (down edge), (11) 
|Mn (e) | ^> 1, localized at y « L y — 1 (up edge). 
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FIG. 2: Energy spectrum of the two dimensional spin-chiral 
ferromagnetic kagome lattice system with fixed boundary un- 
der different spin-chiral parameter (a) <f>=0, (b) <f>=ir/3, and 
(c) <f)=n/2. The shaded areas are the energy bands and the 
lines are the spectrum of the edge states. The solid and dot- 
ted line mean that the edge state is localized near y ~ 1 and 
y ~ L y — 1, respectively. 




FIG. 3: (Color online) Two sheets (Riemann spheres) with 
3 cuts which correspond to the energy bands of the system. 
The RS of the Bloch function is obtained by gluing the two 
spheres along the arrows near the cuts. 



Because the analytical derivation is very difficult, we 
now start a numerical calculation from Eq. @ and draw 
in figure 2 the energy spectrum of the fixed boundary 
system as a function of k x for three values of spin-chiral 
parameter </>. The number of sites A (or B, C) in the y 
direction is chosen to be L y =50. Clearly, one can see that 
the edge states occur in the energy gaps or at the band 
edges. From Fig. 2(b) and (c), one can clearly observe 
that in the spin-chiral cases, i.e., <f> ^ (or 7r), there 
are three dispersed energy bands (the shaded areas) with 
two edge states (the lines) lying in each energy gap. This 
feature is different from that in the case of the square 
lattice in an external magnetic field [TH ] . in which each 
gap has only one edge state. The reason for this difference 
is that the factory in My(e) is no longer a linear function 
of e in the present case. 

Before studying the Hall conductance of this system, 
we simply review the winding number [ill [l3| . which 
is as well as the Chern number a well-defined topological 
quantity. Let us consider the bulk Bloch function at sites 
with y-coordinate of L y . For the Bloch function, and 
'Fq^ compose an eigenvector of M with the eigenvalue p, 



M(e) 



(*>) 



vF 



(>>) 



0) 



(b) 



(12) 



We extend the energy e to a complex energy to discuss 
a wave function of the edge state. Here we use complex 
variable z for the energy. From Eq. (|12p we get 



A(z)- VA 2 (z)-4 



and 



l(*) = - 



Mu(z)+M 22 ( Z )-yg(z) 



(13) 



- 7 M 21 (z), 



-M n (z) + M 22 (z) + VA 2 (z)-4 

(14) 

where A (z) =Tr[M (z)] and $1 = 1 used. Since the 
analytic structure of the wave function is determined by 
VA 2 (z) -4, we consider the RS of a hyperelliptic 



curve uj 2 —A 2 (z) — 4. To make the analytic structure of 
u) to be unique, we have to specify the brunch cuts which 
are given by A 2 (z) — 4 < at Zz=0. Since this condition 
also gives the condition for |p|=l, the branch cuts are 
given by the three energy bands. Therefore, A 2 (z) — 4 
can be factorized as 



u = v/a 2 (z) - 4 = yn i=1 (* ~ A *)' 

where denote energies of the band edges. The RS 
is obtained by gluing the two Riemann spheres at these 
branch cuts along the arrows (see Fig. 3). The Riemann 
spheres are obtained by compactifying the |z|=oo points 
to one point. After the gluing operation, the surface is 
topologically equivalent to the surface shown in Fig. 4. 
In the present model, the genus of the RS is 5=2, which is 
the number of energy gaps. In this way, the wave function 
is defined on the genus-2 RS S g= 2(fca;)- The branch of the 
function is specified as A 2 (z) — 4 > (z — ► — oo on the 
real axis of R + ). Then if z lies in the jth gap from below 
on the real axis (notice that there are two real axes), 
a ^A 2 (z) - 4 > 0, z (real) on R a (a=+, -). So, 

at the energies of the edge states fij , 



V^ 2 (^)-4 = « - MaafoOl (15) 
(/ij GiT, a = +,-). 

In addition, by simple calculation, we can also obtain 



A(e) 



< —2 for j odd 
> 2 for j even 



(16) 



where the energy e (on R^) is in the jth gap. From 
Eqs. (fT5|) . (HHJ) and (fTTj) , we can get that when the zero 
point is on the upper sheet of the RS, the edge state is 
localized at the down edge; when the zero point is on the 
lower sheet of the RS, the edge state is localized at the 
up edge. 

In Fig. 4, on the RS, the energy gaps correspond to 
circles around the holes of the T, g= 2(k x ) and the energy 
bands correspond to closed paths on S g= 2(fc x ). The 
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FIG. 4: (Color online) The RS of the Bloch function. C (fij) is 
a loop formed by the trace of the zero point of &L y -i (z). The 
energy bands are shown by closed loops. The corresponding 
winding numbers are I (C (A i j)) = l for all j. 



Bloch function is defined on this surface. \& 



(b) 



has 



always 2<?=4 zero points at the edge-state energy 
{^L^-i (Mi) = 0)- Since there are two real axes on the 
^g=2{k x ), there are eight fXj's on the RS. However, only 
one of every two gives a zero of 

Changing fc K in one period, we can consider a family 
of S„ = 2(fc a; ). T, g= 2(k x ) can be modified by this change 
yet all the T, g —2(k x ) with different k x 's are topologically 
equivalent if there are stable energy gaps in the two- 
dimensional spectrum. By identifying the topologically 
equivalent Ti g= 2{k x ), we can observe that the fJ-j (k x ) 
moves across the holes and forms an oriented loop C (pj). 
Note that the present case is prominently different from 
the previous case of the square lattice under an external 
magnetic field [ll| . in which fij moves around the holes. 
The two edge-state energy loops C's in the same energy 
gap are tangent at one point and their appearance shows 
a "oo" structure, as shown in Fig. 4. 

As pointed out by Hatsugai [ll|, [l3| , when the Fermi 
energy e F of the two-dimensional system lies in the ith 
energy gap, the Hall conductance is given by the winding 
number of the edge state, which is given by the number 
of intersections I (cti,C (fii)) (= 7(C(^j))) between the 
canonical loop ai on the RS and the trace of fii. In 
the present model, because there are two edge states, 
and, correspondingly, there are two canonical loops in 
one energy gap, the Hall conductance by its definition 
can be written as 



i.edge _ 



^I{C{^j-x)), e F <e Tj 
-^(C(M2i)), e F >e Tj 



(17) 



where exj is the energy at the tangent point in the jth 
energy gap. Similarly, this expression can be obtained 
by the Byers and Yang's formula [20] : Suppose that one 
increases an external magnetic flux $ from to 1 adia- 
batically. According to the Laughlin-Halperin argument 
[2ll I22I ] , when the Fermi energy lies in the jth energy gap 
and e F < e T] [or e F > e Tj ], J(C(^ 2 j-i)) [or I (C (/j, 2 j))] 
states are carried from the down edge (y = 1) to the 
up edge (y = L y — 1) in net. The energy change dur- 
ing the adiabatic process is AE=I (C (/i2j-i)) (— e) V y [or 



AE=I (C (fJ.2j)) (— e) V y ], where V y is a voltage in the y 
direction. This gives the Hall current I x as follows 



AE 



(18) 



where <&o=hc/e is the flux quantum. Then we get an 
expression for o- xy sc as Eq. (fTT|) . 

On the genus g=2 RS, the first homotopy group is 
generated by 4<?=8 generators, ai and f3i, i=l,--- ,4. 
We can observe that (J>2j—i (p-2j) moves one time across 
the jth hole, that means C (/i2j-i) [C (^2j)) ~ Pj and 
I J (C {n2 3 -i)) I [\I (C (p-2j)) |] =1- Considering winding 
direction (See Fig. 5 in Ref. [HI), one can obtain that 
I(C(jJti))=l when € (0,tt), while I (C (m))=-l when 
4> S (— 7T,0) for all i. So, 7(C (/ii))=sgn(sin0), and 



l.cdgc 
'xy 



2, edge _ 
'xy 



— — sgn (sin < 



(19) 



Now we turn back to make an analysis of Fig. 2 with 
the help of the above results. At <p—0, the lower energy 
band becomes dispersionless (Fig. 2(a)), which reflects 
the fact that the 2D kagome lattice is a line graph of 
the honeycomb structure [HI]. This flat band touches 
at k x =0 with the middle band, while the middle band 
touches at k x —^-, ^ with the upper band. So there are 
no bulk energy gaps and the Hall conductance is zero in 
this case. At <f> ^ 0, it the 2D kagome lattice has spin chi- 
rality, and there occur two bulk band gaps, as shown in 
Fig. 2(b) for 0=tt/3 and Fig. 2(c) for 0=tt/2. Then ac- 
cording to Eq. (fH)]) . one can obtain that when the Fermi 
energy lies in the bulk gaps, cr^ dgc =cr^j y cdsc =— Note 
that in the case of </>=§ , the middle energy band becomes 
flat (Fig. 2(c)) due to the particle-hole symmetry. 

Finally let us compare cr° dgc (Eq. (fT9|)) in the 
present model with that in the bulk theory [H, Hoj |. 
In the latter, the bulk Hall conductance has been 



derived to be <r3.;J 5and = 

y 

_3,band e_ 

xy 



^-sgn(sin^), a 



2,band = 
X y 



-0, and 



^-sgn(sin(/>) for the three bands. So, when the 
Fermi energy e F lies in the ith energy gap, the bulk Hall 
conductance cr^ ulk is given by 



i, bulk 

xy 



E j,band 



-sgn(sin ( 



(20) 



where £=1,2. Comparing Eqs. (|20|) with Eq. (|T9[) . one 
can obtain that 



edge _ -bulk 
xy xy 



(21) 



This conclusion keeps consistent with the recently estab- 
lished common recognition on the Hall conductance in 
the systems with and without edges. 

In summary, we have investigated the effect of the 
chiral edge states on the quantum Hall conductance in 
the 2D kagome lattice with edges. According to our de- 
rived Harper equation, there are two edge states lying in 
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each energy gap. They are tangent at one point in the 
gap, thus showing a "oo" structure. The energy loops 
for these two edge states move across the holes in the 
RS. We have also analyzed the winding number of these 
two edge states, which gives the quantum Hall conduc- 
tance <J°y ec —~ ^-sgn(sin</>) when the Fermi energy lies in 



the bulk gap. This conclusion keeps consistent with that 
based on the topological bulk Chcrn-numbcr theory. 
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